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ABSTRACT 

This paper presents galaxy source counts at 24 microns in the six Spitzer Wide-field InfraRed 
Extragalactic (SWIRE) fields. The source counts are compared to counts in other fields, and to model 
predictions that have been updated since the launch of Spitzer. This analysis confirms a very steep 
rise in the Euclidean- normalized differential number counts between 2 mJy and 0.3 mJy. Variations in 
the counts between fields show the effects of sample variance in the flux range 0.5-10 mJy, up to 100% 
larger than Poisson errors. Nonetheless, a "shoulder" in the normalized counts persists at around 3 
mJy. The peak of the normalized counts at 0.3 mJy is higher and narrower than most models predict. 
In the ELAIS Nl field, the 24 micron data are combined with Spitzer-IRAC data and five-band optical 
imaging, and these bandmerged data are fit with photometric redshift templates. Above 1 mJy the 
counts are dominated by galaxies at z < 0.3. By 300 microJy, about 25% are between z ~0.3-0.8, 
and a significant fraction are at z ~1.3-2. At low redshifts the counts are dominated by spirals, and 
starbursts rise in number density to outnumber the spirals' contribution to the counts below 1 mJy. 

Subject headings: infrared: galaxies - galaxies: evolution - star:formation - galaxies: starburst - 
cosmology: observations 



1. INTRODUCTION 

Galaxy counts from IRAS first hint ed at strong evolu- 
tion of luminous infrared galaxies (iHacking et alJ fT987t 
ILonsdale fc Chokshil[l993 iGregorich et al.lll995h . Mid- 
infrared surveys made with the ISOCAM instrument 
at 15 microns revealed an excess population at z~0.8 
(jElbaz et al.lll999l l2002t iFranceschini et al.ll200ll ). This 
population was detected in part by the redshifting of the 
strong 7.7 micron PAH feature into the 15 micron band. 

The high sensitivity a t 24 microns of the MIPS instru- 
ment ([Rieke et al.ll2004T ) onboard Spitzer allows infrared- 
luminous galaxies to be traced to even higher redshifts. 
Sensitivity to obscured starbursts at z~2 is enhanced due 
to redshifting of the 7.7 micron feature into the 24 mi- 
cron band. First source counts have been presented by 
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ICharv et al.l ([2004) for the deep Great O bservatory Ori- 
gins D eep Survey (GOODS) test field, bv lMarleau et al.l 
(pOOl for the First- Look Survey (F L S) an d the GOODS 
test field, and by iPapovich et alJ (|2004l ) for a range 
of Guaranteed Time Observer (GTO) fields plus the 
GOODS test field. These counts show a steep bump in 
the normalized source counts as fluxes decrease from sev- 
eral mJy to about 300 /iJy, that is attributed to the popu- 
lation discovered by ISOCAM at z^0.8 plus a previously 
unseen population of galaxies at z^l-2. At still fainter 
fluxes, the no rmalized counts drop of f to the Spitzer con- 
fusion limit (iRodighiero et all |2006| ICharv et all 12004 
iMarleau et al.ll2004j : IPapovich et al.N2004f) . 

This paper presents a detailed analysis of 24 micron 
source counts from the large fields observed for the 
Spitzer Wide-area InfraRed Extragalactic (SWIRE) pro- 
gram. The SWIRE survey is the largest of the extra- 
galactic Spitzer Legacy Science programs, comprising 
mapping of 49 square degrees at wavelen gth from 3.6 
to 160 microns (|Lonsdale et al.ll2003ll200l . This survey 
is designed to dramatically improve our understanding 
of galaxy evolution, including the history of star forma- 
tion, the assembly of stellar mass in galaxies, the nature 
and impact of accretion processes in active nuclei, and 
the influence of environment on these processes at angu- 
lar scales up to 3 deg (such as the study of clus t ering of 
ultraluminous infrared galaxies bv lFarrah et al.1 (|2006h ). 
The SWIRE dataset encompasses six fields and so is an 
excellent resource for quantifying the effects of sample 
variance. 

The goal of this paper is to present accurate Spitzer 
24/im counts in the flux ranging from 50 mJy to ~300/iJy. 
A detailed model-based interpretation of the counts is 
beyond the scope of this paper. However, the composi- 
tion of the counts is explored for the ELAIS Nl field 
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TABLE 1 

Area and photometric uncertainty for each field 



Field 




iize 


Area 






[deg 


X deg] 


[deg 2 ] 


(microjy) 


ELAIS Nl 


3.1 deg 


X 3.0 deg 


9.64 


38.8 


ELAIS N2 


1.9 deg 


X 2.5 deg 


4.87 


37.6 


Lockman 


3.7 deg 


X 3.0 deg 


11.03 


41.8 


CDFS 


2.6 deg 


X 3.0 deg 


7.67 


38.5 


XMM-LSS 


3.1 deg 


X 3.0 deg 


8.34 


48.2 


ELAIS SI 


2.3 deg 


X 3.0 deg 


7.03 


51.2 



using the IRAC data at 3.6-8^m (|Surace et alJ l2007h 
as well as high-quality optical imaging. This merged 
dat as et provides insights into galaxy SEDs and popula - 
tions (| Rowan- Robinson et "aTll2005t iPolletta et alT 2006) , 
and is used in this paper to break down the 24 /im source 
counts by population and redshift. 

The paper is organized as follows: The SWIRE obser- 
vations and data processing are described in Section [2 
The derivation and validation of source counts are pre- 
sented in Section [3j Comparison with models, including 
recent models updated using Spitzer counts, is in Sec- 
tion |U For the ELAIS Nl field, optical identifications, 
photometric redshifts, and template-fitting lead to source 
counts subdivided by redshift and galaxy type in Section 

El 

2. SWIRE 24 MICRON OBSERVATIONS AND DATA 
PROCESSING 

2.1. Observations 

The Spitzer mapping observations of the six SWIRE 
fields were carried out between December 2003 and De- 
cember 2004. For the mapping at 24, 70 and 160 ^m, 
the fol lowing observing strategy using the MIPS instru- 
ment (jRieke et alJl2004j ) was implemented. The MIPS 
Scan Map Astronomical Observing Template was used 
with medium scan rate. The spacing between scan legs 
was 148 arcseconds, approximately half of the width of 
the 5'x5' field of view (FOV) of the 24/im detector ar- 
ray. Three-degree-long scan legs were used for most As- 
tronomical Observing Requests (AORs). Five scan legs 
were typically used per AOR, resulting in a cross-scan 
size of about 25'. The AORs were arranged into two 
coverage sets, separated by half of an FOV in the cross- 
scan direction. Adjacent AORs within a coverage set 
were overlapped sufficiently to ensure no gaps in cover- 
age of the 70/xm array. The approximate dimensions of 
the maps are given in Table [TJ With this observation 
layout, the typical coverage per point is 44 Basic Cali- 
brated Data (BCD) images, each with exposure time of 
4 seconds, for a total of 160 seconds of integration time 
per point at 24^m, 80 seconds per point at 70/im, and 
16 seconds per point at 160/im. Overlap between rotated 
scans yielded a higher coverage in portions of each map. 

Some regions received different coverage levels than 
described above. A 0.5 square degree region of the 
Lockman field was observed as part of Legacy valida- 
tion operations, for a total integration time about 1.5 
times the nominal ones described above. In the Lock- 
man and CDFS fields, the 24 /im maps include regions 
of higher coverage around the Guaranteed Time Obser- 
vation fields, where scans were initially embargoed. The 



ELAIS SI field received only one coverage, and conse- 
quently half the nominal observing time per point, due 
to observing time constraints. Part of the ELAIS Nl 
field, which was missed in January 2005 due to an obser- 
vatory standby event, was filled in in July 2005, leading 
to some irregularities in field shape and coverage. We 
have accounted for different coverage levels and sensitiv- 
ities in our analysis. 

2.2. Data processing 

The 24 micron raw data were processed by the S10.5 
version of the SSC pipelines. The flat-field images used 
were made separately for each scan mirror position, us- 
ing SWIRE scan map data. The SWIRE processing be- 
gan with the BCD (Basic Calibrated Data) images. For 
all fields, to even out variations in the background, the 
median of all the pixels in each image was subtracted 
from each pixel in that image. In the XMM-LSS field, 
removal of dark-latent artifacts was necessary, owing to 
the proximity of the very bright source Mira to the MIPS 
mapping region. For each scan leg, a self-calibration 
flat image was produced from the median of all uncor- 
rected BCDs in that leg, and was normalized and divided 
into those BCDs before the background-subtraction step. 
This technique was successful because the intensity of 
these latents did not change appreciably over the 30- 
minute duration of each scan leg. 

To make maps, the corrected BCDs were coadded 
into large mosaics, using the SSC's MOPEX software 
(jMakovoz fc Marleadl2005T ). The mosaics were all made 
with 1.2" pixels. Cosmic rays were removed using multi- 
frame outlier rejection with a threshold of 2.5 a. 

Source ex traction was performed on the mosaic using 
SExtractor (|Bertin &: Arnoutlll996l ). A local background 
of size 128x128 pixels was computed for the maps, and all 
the noise calculations were weighted by the inverse square 
of the coverage map. Photometry was output as mea- 
sured in several apertures and by various extended source 
techniques. For the present study, we have used photom- 
etry in a 5. 25" radius a perture for point sources, and the 
Kron flux (jKronl H~980| ) for extended sources. Aperture 
corrections were computed from comparing smaller aper- 
tures to a 30.6-arcsecond-diameter aperture. Then, an 
additional correction of 1.15 was applied to both aper- 
ture and Kron fluxes, to match the procedure used by 
the SSC pipelines. This additional factor is derived from 
the fraction of light outside the 30.6" aperture in theoret- 
ical TinyTim models of the point spread function. The 
resulting overall correction for the 5.25" aperture is 1.78. 

The calibration supplied with the data as- 
sumes a source spectrum o f a 10,000 K blackbody 
(jMIPS Data Handbook! |2006[ K Galaxies detected at 24 
lim generally have an S EP slope of F v ~ v^ 1 or redder 
(R owan-Robinson et al.l l2005f ). Accordingly, we have 
color-corrected our photometry to an assumed source 
spectrum of F v ~ v^ 1 by dividing by a factor of 0.961. 
Additionally, we have multiplied by a factor of 1.018 
to account for a small change in the 24 /xm calibration 
facto r made in December 2005 (jMIPS Data Handbook 
120061) . 

For this work, sources are assigned the Kron flux if 
the SExtractor ISO-AREA value is at least 100 pixels, 
and the stellarity is less than 0.8. All other sources are 
considered point sources, and the 5.25" radius aperture 
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measurement is used. 

2.3. Data quality assessment 

Estimating photometric uncertainties for extractions 
from mosaicked images is complicated by correlations be- 
tween pixels. To estimate the uncertainty in each field, 
we measured fluxes in thousands of randomly placed 
apertures that fell in regions of typical coverage. Gaus- 
sians were fit to the distribution of fluxes in these aper- 
tures, yielding robust uncertainty estimates. The pho- 
tometric uncertainties in the 5.25" aperture are given in 
Table [TJ The areas in each field used in this work are 
also given. The relatively high noise value for the XMM- 
LSS field is due to its low ecliptic latitude, leading to a 
high zodiacal background relative to the other SWIRE 
fields. The noise in ELAIS SI is higher because this field 
received half the integration time of the other fields. 

We have estimated completeness by inserting simu- 
lated sources into the ELAIS Nl, ELAIS N2 and XMM- 
LSS fields, and computing the fraction that are recov- 
ered. The completeness curve shifts to higher fluxes for 
XMM-LSS, approximately following the ratio of noise 
levels in Table [TJ We have accordingly shifted the ELAIS 
Nl curve by the ratio of noise levels to estimate complete- 
ness in Lockman, ELAIS SI and CDFS (Figure [TJ). In 
the best fields, 70% completeness is attained for flux den- 
sities above 0.35 mJy. Based on the noise estimates in 
Table [TJ this flux density corresponds to a signal-to- noise 
ratio of about 7. Ordinarily the completeness would be 
expected to be higher at this SNR, but in our case the 
SExtractor detection parameters used have limited the 
completeness more than the noise has. 

A direct estimate of completeness may be made by 
comparing the SWIRE data with deeper data. We have 
made such a comparison by processing GTO data in the 
CDFS field using the same methods used for the SWIRE 
data, and matching the extractions with a 3.6" search 
radius. The completeness obtained from this comparison 
is in good agreement with the simulation results. 

The comparison to deeper data is also useful as an 
assessment of reliability. The ratio of all matched sources 
to all extracted SWIRE sources is also plotted in Figure 
[TJ The reliability defined in this way is better than 98% 
in all bins down to 0.35 mJy. 

3. SWIRE NUMBER COUNTS 

Although the contribution of stars to the 24 mi- 
cron number counts is small at faint (sub -mJy) fluxes 
(jMarleau et al J 12004 : [Fapovich et all 120041) , it must be 
taken into account at fluxes of several mJy and brighter. 
We have used 2MASS data to identify stars. Our 24/im 
extractions were matched against the 2MASS catalog us- 
ing a 3" search radius. The resulting K s — [24] color- 
magnitude diagram is shown in Figure [2] [24] is the 
24fim magnitude, where we have used a zeropoint of 
7.43 Jy in our color-corrected flux units (correspond- 
ing to 7.14 Jy for a 10 , 000 K blackbody spectrum, see 
IMIPS Data Handbook! ([20061 )). The 3a limit for our 
most sensitive fields corresponds to [24] = 12.0, well 
above the 2 MASS completene ss limit in unconfuscd sky 
of K s ~ 14.3 ([Skrutskie et al.ll2006f ). For our high Galac- 
tic latitude fields, the 2MASS survey is sensitive enough 
to detect all stars in the SWIRE 24 /im catalogs. The 
dashed line in the figure shows the simple color cut used 



TABLE 2 
Total SWIRE 24 fim counts 



Average F v 


Normalized counts 


Uncertainty 


(mJy) 


(mJy ° deg ) 


(mJy ° deg J 


0.405 


1018.81 


4.93 


0.510 


912.33 


5.34 


0.642 


817.73 


6.15 


0.809 


722.46 


6.82 


1.018 


645.91 


7.81 


1.282 


592.51 


8.97 


1.613 


509.69 


9.08 


2.279 


516.15 


10.21 


3.612 


547.15 


13.38 


6.423 


484.33 


17.76 


12.816 


402.43 


26.17 


28.691 


364.24 


41.77 


72.069 


299.29 


64.37 



to identify stars for the purpose of this st atistical study. 
A sim ilar color discriminant was used by IVaccari et al.l 
([2005] ) with ISO 15 pirn fluxes. 

Integral counts for the stars and galaxies are shown in 
Figure [3J The integral counts here are the total in each 
field, without normalization by the areas in Table[TJ The 
contribution of stars and galaxies to the number counts 
is approximately equal at 30 mJy. 

Euclidean-normalized number counts were obtained by 
computing differential counts with each source individ- 
ually weighted by a factor of flux raised to the power 
2.5. We have plotted these counts for all fields in Figure 
[4j A completeness correction has been applied for each 
field, using the curves shown in Figure [TJ shifted for each 
source by the square root of the ratio of the coverage to 
the mean coverage. Where the completeness correction 
becomes large, as in the faintest bin shown for ELAIS 
SI, the counts are under-corrected. 

The field-to-field comparison shows differences on the 
order of 10% in the submJy region. At these fluxes there 
are of order 10,000 galaxies per field (Figure [3]) so that 
counting statistics can account for only a fraction of this 
difference. Even more striking is the large variations in 
the vicinity of 3 mJy. Taken together, these differences 
are almost certainly due to sample variance, significant 
even on 3-degree scales. 

Total SWIRE counts were computed from the com- 
bined sample with corrections for incompleteness ap- 
plied. Uncertainties in the counts were estimated using 
a bootstrap technique, in which samples of identical size 
to the original are drawn with replacement from the to- 
tal sample. The standard deviation of the counts derived 
from the bootstrap samples provides the error bars. The 
uncertainties are generally quite similar to what is ex- 
pected from Poisson statistics, but are somewhat larger 
at higher flux levels. Since the counts are derived from 
six widely spaced fields, sample variance effects should 
be averaged out. The total SWIRE counts are shown 
in Figure [5] and are tabulated in Table [2j We have also 
plotted the c ounts for the FLS dMarleau et al.ll2004l ) and 
GTO fields (|Papovich et~aTl l2004h for comparison. To 
put the Papovich counts on the SWIRE calibration, we 
multiplied the fluxes from that work by a factor 1.059 
to account for our color correction, the change in cali- 
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^(24/xm) [mJy] 

Fig. 1. — Filled symbols: Differential completeness from simulations for ELAIS Nl, ELAIS N2 and XMM-LSS fields. The curves for the 
other three fields were obtained by shifting the EN1 curve as described in the text. Open squares: Differential reliability defined as the 
fraction of SWIRE sources found in the deeper GTO data in CDFS. 




Fig. 2. — Color-magnitude diagram of 2MASS K s and [24] for 
all six SWIRE fields. The dashed line shows the color cut used to 
separate out the stellar contribution (below the line) to the 24/im 
number counts. 



bration factor (|MIPS Data Handboolj|2006D , and a very 
slight difference in aperture corrections. The FLS counts 
were shifted assuming a flux scale fainter than ours by 
a factor of 1.11, due to the color correction, the change 
in calibration factor, and our 5% larger aperture correc- 
tion for a 30.5" diameter aperture. The steep rise in the 
normalized counts found by these works is confirmed by 
our results. There are some differences at brighter flux 
levels, but these appear to be within the range of the 
SWIRE field-to-field variations shown by the error bars. 
Also of note is a "shoulder" at about 3 mJy, seen in the 
total SWIRE and lPapovich et all (|2004D counts, though 



not in the Ma rleau et a l. (2004) counts, and also not in 
the ELAIS SI counts (Figure [J). 

4. COMPARISON OF COUNTS WITH MODELS 

A plot of the normalized counts with several mod- 
els overlaid is in Figure El The count s used include 
the SWIRE counts, the iPapovich et~aTI (|2004D counts 
at fainter fluxe s , and IRAS counts from the sample 
of IShupe et all (|1998l ) corrected to the mean MIPS 
wavelength of 23.7 /im. The fitting of several pre- 
.Spiteer-launch mo dels have already been di scussed in 
ICharv et all (|2004l ) and lPapovich et alJ (|2004h . Here we 
focus on more recently availa ble models. One of these 
models is described in detail in lLagache et al.l (|2004f ) and 
is not outlined here. We give a thumbnail sketch of the 
other models in turn. 

4.1. Details on individual models 



The iGruppioni et alJ (120051) model i s base d on star- 
burst +cirrus gala xies from Pozzi et al.l (|2004f ) and AGN 
components from lMatute et aTp2006f)~ These works are 
15 micron models based on data from the ELAIS sur- 
vey (jRowan- Robinson et al.ll2004h . They have been con- 
verted to 24 microns using the appropriate ratio of 24/im 
to 15/xm of each population. Though this work was pub- 
lished after the launch of Spitzer, it did not involve an a 
posteriori correction. 

The predicted 24/im counts by iFranceschini et aTl 
(2007) are based on an IR multi-wavelength evolution- 
ary model fitting a variety of data from Spitzer, ISO, 
SCUBA, and reproducing data on the diffuse far-IR back- 
ground observed by COBE. The model basically adopts, 
in addition to normal spiral galaxies and type-1 quasars, 
two populations of starbursting galaxies with different lu- 
minosity functions and evolution rates. One population 
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F„(24/im) [mJy] 



Fig. 3. — Integral counts for the contiguous regions in each SWIRE field. Galaxy counts are the symbols connected by lines, while stellar 
counts are unconnected. Note that the counts have not been normalized by area so as to provide some separation between the distributions. 



consists of moderate-luminosity star/bursts with peak ac- 
tivity at z~l (the LIRG population), the second contains 
more luminous objects (the ULIRG and HYLIRG) with 
maximum activity phase at z^2. This evolutionary bi- 
modality is apparently required for a combined fit of the 
multi- wavelengt h coun ts and redshift distributions. 

The IXu et all (|2007f ) model includes a new evolution 
model for dusty galaxies a nd the E2 model for the E/SOs 
evolution (|Xu et alj 120031 ).. The new dusty galaxy evo- 
lution model is a modified version o f mod el SI, the 
starburst-dominant model, in lXu et" al. (2003). The def- 
initions of galaxy populations (starbursts, normal-disks 
and AGNs) are the same as in model SI, as are the evo- 
lution parameters of the normal-disks and AGNs. Both 
the luminosity evolution rate and density evolution rate 
of the starburst galaxies are on the order of (1 + z) 3 for 
z < 1 and (1 + z)~ 15 for larger redshifts. The major 
changes are in the SEDs of the starbursts and normal- 
disks: For both populations the strength of the PAH 
features in the wavelength range of 5/zm < A < 12/zm is 
raised by a facto r of 2. 

The model of iRowan- Robinson ( 20 0?]) is a modifica- 
tion of the counts model of IRowan- Robinson! (|2001l ). It 
retains the four infrared SED types of the latter, with 
no modification of th e SE Ds (unlike the n ew models of 
iLagache et all (|2004l ) and IXu et all (|2007l )). The four 
types still all undergo strong luminosity evolution but 
the evolution is allowed to be at a different rate in the 
different components. The evolution of the cirrus (qui- 
escent) component is now shallower than the evolution 



of the M82 starburst component, and is considerably 
steeper for the Arp 220 (high optical depth starburst) 
component. The evolution function is also modified so 
that the evolu t ion is less steep at z < 0.5 (as in the 
ILagache et al.l (|2004l ) model). The stronger evolution 
in the Arp 220 component make s the mode l simil ar to 
those of iFranceschini etldl (Il997f) . IXu et al.l (|1998f ) and 
iDole. Lagache fc Pugetl 1)20031 ) . but it still has the dis- 
tinct feature of strong evolution even in the quiescent 
component. 

4.2. Comments on model fits 

It is difficult to fit a steep bump in the normalized 
counts without restricting the luminosities in the models 
to a narrow range. That said, the models are successful 
in peaking at about the same flux level (0.3 mJy) as the 
source counts. The counts rise rather more steeply with 
decreasing flux, however, so that they peak with counts 
about 25% higher than the models predict. Some of this 
difference is due to the color correction and calibration 
changes applied to t he counts, which inc rease flux by 
about 5.9% for the (|Papovich et alj l2004h counts, and 
increase normalized counts by about 13%. The model 
curves could be scaled up by this factor, but in many 
cases would overpredict brighter sources. 

Models which allow evolution of only the more lu- 
minous sources, including the possibility of discontinu- 
ous changes in the evolution rate, or large changes to 
source SEDs, have less difficulty in fitting the counts. 
By contrast, models which use smoothly-changing evo- 
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Fig. 4. — Normalized differential number counts for the six SWIRE fields. For clarity, representative Poisson error bars are shown at the 
top as computed for the ELAIS Nl field. 



lution rates, like those of I Rowan- Robinson! (|2001f ) and 
iRowan-Robinsonl (|2007l ), have more difficulty. It is im- 
portant to establish what really is needed to account for 
the observed counts. It may be that PAH emission an d 
silicate absorption features at z ~ 2 (Houck ct aLl[2005TI . 
or even silicate emission from AGN (jHao et all 120051 ) , 
are stronger or play a larger role than is assumed by the 
models. 

5. COMPOSITION OF 24/iM COUNTS IN ELAIS Nl 

We now turn to an investigation of the composition of 
the 24/tm counts. Here we rely on the full optical and 
Spitzer-IR AC data available in the ELAIS Nl field. 

5.1. Optical identifications, photometric redshifts and 
infrared template fits 

We have associated our band-merged Nl catalogue, 
which contains all 24 /zm detections with a 3.6/im 
counterpart to ensure high reliability, with the optical 
UgriZ catalog generat ed from the INT WFS data by 
IBabbedge et"alT (|2006) , using a search radius of 1.5 arc- 
sec. This search radi us is appropriate for SWIRE s ources 
detected at 3.6 /mi ()Rowan-Robinson et al. 2005). The 
proportion of 24 /im blank fields brighter than S(24) = 
200 /iJy, to the WFS survey limit of r - 23.5, is 28%. 

All sources with optical associations have been 

run through t h e ph otome tric redshift c ode of 

IRowan-Robinsonl (|2003f ) and IBabbedge etaD (|2004f ). 
with small modifications to the spectral energy distri- 
butions (SEDs) described by iRowan-R .obinso n et al.l 



(2005). 20% of 24 /im sources brighter than 200 /tJy 
with optical associations failed to get a photometric 
redshift because there are less than 4 bands detected 
at UgriZ, 3.6, 4.5 /tm or because the reduced x 2 1S too 
poor (> 10). This can be because of erroneous optical 
photometry, incorrect optical associations, or because 
the range of optical templates used fails to characterize 
all sources. The latter possibility is worth further study, 
especially since, for example, strongly reddened quasars 
or galaxy+quasar combined SEDs are not included in 
the SED library. 

Finally, the SEDs of all 24 /im sources with pho- 
tometric redshifts are fitted with infrared templates 
added to the optical/near IR galaxy and quasar tem- 
plate used in the photometric redshift fit, provided 
there is an excess in at least two bands relative to 
the optic al /NIR galaxy model. The te mplates used are 
those of iRowan-Robinson et al.l (|2004f ). also used ear- 
lier in IRowan-Robinsonl (|2001h : cirrus, Arp 220 star- 
burst and a mixture of M82 starburst and AGN dust 
torus. E xamples of SED fits using thes e templates are 
given by IRowan-Robinson et al.l (|2005l ). For sources 
with an infrared excess only at 24/zm, for which we 
can not at this stage characterize the far infrared SED, 
we have fitted an M82 starburst template. Photomet- 
ric redshift distributions of 2 4/tm sources are shown in 
IRowan-Robinson et all (|2005l ). 

5.2. 24 fJ-m counts subdivided by redshift and infrared 
template type 
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Fig. 5. — Plot of normalized number counts from the average of all SWIRE fields, along with those of IPapovich et al.l ll200l) and 
IMarleau et~aTI II2004I) . 



The results of the template fitting may be used to break 
down the infrared counts by redshift range and type. The 
fraction of counts per redshift range is shown in Figure 
El A r elated plot is shown in Figure 2 of [Babbcdg e et aD 
(2006) but here we show the results as a fraction of to- 
tal counts. Overlaid are the same fra ctions taken from 
the m odels of lLagache et al.l (|2004D and lRowan- Robinson! 
(|2007f ). The templates fits and the model are generally in 
agreement for z < 0.3. The model yields a higher fraction 
in the higher redshift bins than the template fits; however 
since some ELAIS Nl sources did not receive a redshift 
in the fitting procedure, this may explain the differ ence. 
Similar results were foundbv lLeFloc'h et al.l (|2005ft and 
iPerez-Gonzalez et al.l (|2005|) in which the models gener- 
ally underpredicted the sources in the low-redshift bins 
and overpredicted the number at higher redshifts. 

Figure [8] shows the 24 /iin differential counts subdi- 
vided by infrared template type. We see that the steep 
rise between 1 mJy and 300 /iJy is caused primarily by 
starbursts, and sources with a infrared excess at 24 fun 
only. To really understand the 24 /jm counts, we need 
70 /im data for the sources with infrared excess over the 
fitted template at 24 /mi only. The fraction of sources 
with quasar-like infrared SEDs is below about 30% at 
bright fluxes, decreasing to less than 15% at relatively 
faint fluxes. 

6. SUMMARY 

The variation between fields of the SWIRE 24/xm num- 
ber counts shows the importance of sample variance. For 



fluxes near a few mJy, the field-to-field variation in the 
counts is 50% to 100% larger than is accountable solely 
from Poisson uncertainties. Based on redshift estimates 
from template fitting and from models, the counts vari- 
ations near 3 mJy may result from large-scale structures 
around z ~ 0.5. At fainter flux levels, the SWIRE 24/mi 
catalogs will be useful for tracing structures at z ~ 1 and 
beyond. 

The 24 /im counts seem exceptionally interesting 
and most models published before the launch of 
Spitzer miss badly. There may be new populations- 
iRowan-Robinson et al.l (|2005T ) claim there is a large pop- 
ulation of Type 2 AGN that have not been accounted for 
in previous models. The SWIRE counts confirm a steep 
rise in the counts from about 2 to 0.3 mJy that is quite 
difficult for models to fit. 

Template fitting to the SWIRE EN1 sources is suc- 
cessful for at least 85% of the galaxies. The photometric 
redshifts indicate that at least a quarter of the sources in 
the 1-10 mJy interval have photometric redshifts in the 
range 0.3 to 0.8. 

The 25/im band onboard IRAS is closest in wavelength 
space to the MIPS 24/tm p assband. IRAS co unts are 
complete to about 250 mJy (|Shupe et al.lll998t ) and the 
SWIRE counts go to about lOOOx fainter. SWIRE re- 
sults should be the best 24 /im counts above 250 /j,Jy 
because it will cover the most area of the Spitzer extra- 
galactic surveys. The full area of Spitzer surveys solidify 
the counts up to about 30 mJy, almost allowing link- 
ing up with the all-sky IRAS counts, but leaving a gap 
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Fig. 6. — Plot of normalized number counts with several models overlaid. 



of a factor of ten in flux. It will be left for the Akari 
(ASTRO-F) and WISE missions to bridge this gap. 

We are grateful to the referee for comments, which have 
led to substantial improvement of this paper. This work 
is based on observations made with the Spitzer Space 



Telescope, which is operated by the Jet Propulsion Lab- 
oratory, California Institute of Technology under a con- 
tract with NASA. Support for this work was provided by 
NASA through an award issued by JPL/Caltech. SB was 
supported by a Fondazione Ing. Aldo Gini 2006 grant. 
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